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Abstract 

Boundary form factor axioms are derived for the matrix elements of local boundary 
operators in integrable 1 + 1 dimensional boundary quantum field theories using the 
analyticity properties of correlators via the boundary reduction formula. Minimal 
solutions are determined for the integrable boundary perturbations of the free boson, 
free fermion (Ising), Lee- Yang and sinh-Gordon models and the two point functions 
calculated from them are checked against the exact solutions in the free cases and 
against the conformal data in the ultraviolet limit for the Lee- Yang model. In the 
case of the free boson/fermion the dimension of the solution space of the boundary 
form factor equation is shown to match the number of independent local operators. 
We obtain excellent agreement which proves not only the correctness of the solutions 
but also confirms the form factor axioms. 



1 Introduction 

The bootstrap program aims to classify and explicitly solve 1+1 dimensional integrable 
quantum field theories by constructing all of their Wightman functions. The first stage is 
the S-matrix bootstrap: the scattering matrix, connecting asymptotic in and out states, 
is determined from its properties such as factorizability, unitarity, crossing symmetry and 
the Yang-Baxter equation (YBE) supplemented by the maximal analyticity assumption. 
The result is the complete on-shell solution of the theory, i.e. the spectrum of excitations 
and their scattering amplitudes, which can be related to some independent definition of 
the model as a perturbed conformal field theory or a Lagrangian QFT (for reviews see 
[HE]). The second step is the form factor bootstrap, which allows one to determine 
matrix elements of local operators between asymptotic states using their analytic properties 
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originating from the already known S-matrix. Supposing maximal analyticity leads to a set 
of solutions each of which corresponds to a local operator of the theory. The form factors are 
then used to build the correlation (Wightman) functions via their spectral representations, 
yielding a complete off-shell description of the theory (see 0111 for reviews). 

The first step of an analogous bootstrap program for 1+1 dimensional integrable bound- 
ary quantum field theories, the boundary R-matrix bootstrap, has been developed for sev- 
eral theories. In boundary theories the asymptotic states are connected by the reflection R- 
matrix, which satisfles unitarity and boundary crossing unitarity; for integrable boundary 
QFT, it also satisfies the boundary YBE (BYBE) and boundary bootstrap requirements. 
These equations supplemented by maximal analyticity assumptions make possible to de- 
termine the reflection matrices and provide the complete information about the theory on 
the mass shell [H]. 

For the second step matrix elements of local operators between asymptotic states have 
to be computed. In a boundary quantum held theory there are two types of operators, the 
bulk and the boundary operators, where their names indicate their localization point. Due 
to the broken translational invariance one point functions of bulk operators may acquire 
nontrivial space dependence behaving analogously to the two point functions in a bulk 
theory. Indeed this one point function can be calculated in the crossed channel, where 
the role of time and space is changed and the spatial boundary appears as a temporal 
one represented as an initial (boundary) state in the matrix element. Inserting a complete 
system of the bulk Hilbert space a spectral representation for the one point functions can 
be obtained in terms of the bulk form factors and the matrix element of the boundary state 
PlIZI. Truncating this expansion at flnite intermediate states provides a convergent large 
distance expansion. However, matrix elements of boundary operators cannot be computed 
in this way and the purpose of the present paper is to develop a technique to compute their 
correlation functions. 

In this paper we initiate the second step of the boundary bootstrap program, namely 
the boundary form factor program for calculating the matrix elements of local boundary 
operators between asymptotic states. We derive their analytic structure from that of the 
R-matrix which, when supplemented by the assumption of maximal analyticity, leads to 
their determination. In the bulk case, it was shown in [H| that the solution space of the 
form factor equations can be brought into one-to-one correspondence with the operator 
content of the model. Based on this, we expect that the classification of the solutions of 
the boundary form factor axioms provides information on the boundary operator content 
of the theory, which in the ultraviolet limit is in a one-to-one correspondence with the 
Hilbert space of the model. Using the explicit form of the boundary form factors the 
spectral representation for the boundary correlation functions can be obtained. 

The paper is organized as follows: first we define the boundary form factors by intro- 
ducing asymptotic in and out multi-particle states, which are related by the multi-particle 
refiection matrix. Simple crossing relations are presented from which the form factor ax- 
ioms follow easily, and then the axioms are verified by some consistency requirements. We 
outline a general strategy to solve theories with diagonal bulk scattering and boundary 
refiection amplitudes, and to compare the resulting two-point functions with their ultra- 
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violet limits. This idea is applied to integrable boundary perturbations of several models, 
such as the free boson model, free fermion (alias Ising) field theory, the scaling Lee- Yang 
model and sinh-Gordon theory. Appendix A contains a heuristic derivation of the cross- 
ing relations from the boundary reduction formula while in Appendix B we present a 
formal derivation of the boundary form factor axioms from the boundary version of the 
Faddeev-Zamolodchikov algebra. 

2 Boundary form factors 
2.1 Definitions 

The Hilbert space of a boundary quantum field theory consists of multi-particle states, 
which can be labeled by the particle species and the corresponding particle energies. To 
simplify the notations we restrict ourselves to theories containing only one particle type 
with a given mass m. In 1+1 dimensions it is convenient to work with the rapidity variable 
9i; the energy Ei of the particle can be written as Ei = m cosh 9i, while the momentum 
is Pi = msinh^j. Following the evolution of the multi-particle state in time to t — — oo 
the particles get far away form each other and from the boundary, therefore forming an in 
state which is equivalent to a free multi-particle state and is denoted as^ 

\9l, 02, ■ ■ ■ , On)in ; 9i> 92> ■ ■ ■ > 9n> Q 

Positivity of all incoming rapidities is a consequence of the assumption that the boundary 
is at the right end of the half line and it is a major difference from the bulk situation. This 
difference is essential because it influences the analyticity domain of matrix elements. 

For t +00 all the scatterings and reflections are terminated, the particles are again 
far away from each other and from the boundary forming the out state, 

l^n • • • 5 (^m)out ; 9]^ < 92 < ■ ■ ■ < 9^ <Q 

which is again equivalent to a free state. By the standard assumption of asymptotic 
completeness, the two sets of states form a complete basis separately and are connected 
by the reflection matrix, which is the boundary analogue of the S matrix. In an integrable 
theory, due to the infinite number of conserved charges, there is no particle creation {n = 
m), the set of rapidities changes only sign 9^ = —9i, and the reflection matrix factorizes 
into the product of pairwise bulk scatterings and individual reflections 

\9l, 92, ... , 9n)in = Y\_^(^i " + ~ ^1' ~^2, • • • , —9n)out (1) 

i<j i 

where S{9i — 9j) connects the two particle asymptotic in and out states in the bulk theory 
(without the boundary) 

^In general, particles in an interacting two dimensional quantum field theory have an efi^ective fermionic 
statistics with the sole exception of free bosonic theories, for which it is necessary to allow equality in the 
ordering of the particle rapidities. 
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I^^i, ^2)^' = 5(^1 - ^^2)1^2, ^i)^*' depicted as 

01 

It is defined originally for 6i > 62 but can be analytically continued for complex rapidity 
parameters such that the extended function (denoted the same way) is meromorphic and 
satisfies unitarity and crossing symmetry 

s{e)s{-e) = 1 , s{t7i -e) = s{e) 

It might have poles on the imaginary axis at locations 6 = iuj with residue —iresg=iu^ S (9) = 
T'j, some of which correspond to bound states. 

The amplitude R{9) connects the one particle asymptotic states in the boundary theory 

-e 

\9)in = R{9) \ - 9) out depicted as 

e 

It can also be extended from the fundamental domain 6' > to a meromorphic function 
on the whole complex 9 plane satisfying unitarity and boundary crossing unitarity 

R{9)R{-9) = 1 , Riin - 9)S{29) = R{9) 

R{9) may have poles at imaginary locations 9 = ivj (0 < Vj < 7r/2), some corresponding 
to excited boundary states. If the interpolating field has a nontrivial vacuum expectation 
value then generally there is also a pole at 6^ = z7r/2 with residue 

-i Res R{9) = — (2) 

The boundary form factor is defined as the matrix element of some local boundary 
operator, 0{t), between asymptotic states 

out {9 1 , 9„...,9JO{t)\9,,92, 

These form factors are defined only for 6'i > 6*2 > • ■ ■ > > and 9[ < 92 < ■ ■ ■ < 9'^ < 0. 
We can introduce other form factors as 

out{9i, ^'2, • • • , 9^\0{t)\ - 9i, -92, . . . , -9n)out = ^2' • • • ' ^m! ~^2, • • • , -6'n) 

and consider them as a continuation of the original ones in the rapidities. Expressing these 
form factors (via the boundary reduction formula IHj) in terms of correlation functions 





• • • ) "n/in 

FZJ9' 9'..., 9- 9,, 92,..., 9n)e~'^'^^ cosh e.-E cosh.;) 



4 



an analytic continuation can be performed for any (even) complex values of the rapidity 
parameters. As a result the generalized form factors are meromorphic functions of the 
rapidity parameters, and we shall assume that their poles always have physical origins 
(maximal analyticity assumption). From the crossing formula 

(^2, ■■■,e^;e^+ i-K, 61,62,..., 6n) + disc (3) 

derived in Appendix A, we can express all the form factors in terms of the elementary form 
factors 

OMi (0 1 C (0)1 6^1, 62,..., 6n)in = F^{6l, 62,..., 6n) 

It is important to notice that the boundary form factors F^{6i, . . . , 6n), in contrast to the 
bulk case, do depend in general on all the rapidities 6i, not just on their differences, since 
in the presence of a boundary Lorentz invariance is broken. 



2.2 Axioms 

In the Appendices we derive all the following properties of the matrix elements of local 
boundary operators valid in any integrable boundary quantum field theory. Following the 
general philosophy in the bulk case [3] we take them as axioms defining the local operators 
via their matrix elements. 
I. Permutation: 



-^jf ■ ■ - i^i) ^i+i7 . . . i6n) = S{6i — 6i^i)F^{6i, . . . , 6'j+i, 6i,..., 6'„) 



Oi 



II. Refiection: 





, 6n~l, 6n) = R{6n)F^{6i, . . . , 6n-l, —6n 





III. Crossing reflection: 

F^iei, 62,..., en) = R{tTT - (2z7r - ^1, ^2, ... , On) 




The singularity structure of the form factors is determined on physical grounds and 
axiomatized as follows: 
IV. Kinematical singularity 




equivalently described as 

Res F^^,{-9 + m, 9',9,,..., 9^) = [r{9) -\[S{9 - 9i)R{9)S{9 + 9^) ) F^{9u ..-A 




V. Boundary kinematical singularity 

-i Res FS^,{e + |, e„ . . . , e„) = - f[ - 9;))f^{9,, ...,9n) 

1=1 



-i Res e 
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VI. Bulk dynamical singularity 



-t ResFf+2(^ + e' -tu,ei,..., On) = TF^^,{e, 01,..., 



G+ii) 



-i Res 




e+i-o 




VII. Boundary dynamical singularity 



-I Resi^f+i(^^i,...,e, 



d=iv 



-i Res 





We note that equations similar to some of ours have been obtained earlier studying 
boundary form factors in specific spin chains. Using a concrete realization for the Hilbert 
space and the operators, these equations were extracted originally for the XXZ and XYZ 
models in [l^ and extended for other spin chains in [TT]. By extending the bulk free field 
representation for the boundary sine-Gordon model the analogues of XXZ equations were 
obtained in [12]. In all these approaches, however, there is no analogue of the axiom V, 
without which the equations do not determine completely the form factors as can be seen 
on the example of the sinh-Gordon model. In contrast, in our approach the form factor 
axioms are firmly established from first principles of local quantum field theory, thus they 
are valid in a general setting. As a further result of our systematic approach the axioms 
found form a complete system ready to be solved. 



2.3 Consistency checks 

Before proceeding to concrete examples we perform a few consistency checks of the axioms. 
First we note that they are self-consistent in the sense that for specific rapidities the n + 2 
particle form factor can be connected to the n particle form factor either by the kinematical 
singularity equations or by using twice the boundary kinematical equations, and the two 
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procedures give the same result. Indeed taking double residue in the first case, first at 
6 = and then at = gives 




Taking now the residue at 6* = ^ first then at 6*' = y and using that S'(O) = — 1 gives 



i Res i Res F^.^i-O + iir, 61', 6*1, ... , On) = 

Q ITT q' ITT 

The two different orders of taking the residues differ by a factor of 2 since in the first case 
after taking the residue at 6* = 9 we get a factor f{26 — in) which has a zero at 9 = 
(due to 5(0) = —1 the bulk minimal form factor vanishes at the threshold: /(O) = 0). In 
the second case after taking the first residue a factor f{9 — z|) appears. When expanding 
around 6* = 2| to take the second residue there appears a factor 2 due to the difference 
in the arguments of this particular factor (all other terms are identical in the two cases). 
Combining the crossing symmetry of the S-matrix with the definition of g ^ the two 
expressions are easily seen to be equivalent. 

It is worth emphasizing that in the boundary kinematical singularity axiom it is the 
particle-boundary coupling constant g which appears although the residue of the reflection 
factor determines only g^. There are known examples where in two physically different 
situations the fundamental reflection amplitudes are the same and the two cases are dis- 
tinguished only by the sign of g (e.g. the boundary Lee- Yang model with 1 boundary and 
with $ boundary with a particular value of the boundary coupling ^Hl ^ see in more details 
in Sec. 3.3). Because of axiom V the solutions of the form factor axioms are different for 
the two cases, as shown in detail in Sec. 3.3. 

As a second test we relate the two disconnected physical domains {in/out) of the def- 
inition of the form factor. By permuting successively each rapidity to the last position, 
applying a reflection and permuting back to their original position we obtain that 

F^{9,, . . . , ^„) = J] 3(9, + 9,) H R{9,) J] 3(9, - 9^)F^{-9,, . . . , -9^,) 

i<j i i<j 

The product appearing is nothing but the multi-particle R-matrix, (|T}, which connects the 
in and out states. 

Finally we use the fact that the reflection matrix can be considered as a special form 
factor (of the identity operator Id) whose analytic properties are well known. By deflnition 

F2"(^' + zn, 9) = out{9'\Id\9)in = R{9) out{9'\Id\ - 9)out = R{9)5{9 + 9') 
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Now using the permutation property and tricks as above we have 

Fi\e' + in, e) = s{in + e' - e)F^'^{e, e' + m) = siin + e' - e)R{n^ + e')5{e + e') 

which, due to the boundary crossing unitarity, is equivalent to the previous expression. 

2.4 General solution 

In this section we describe the general procedure we use to obtain the solutions of the form 
factor equations in the various specific models. In doing so we emphasize the similarities 
and the differences between the boundary and bulk form factors and also separate the 
(boundary) operator dependent parts from the ones that depend on the specific field theory 
considered but are independent of the operators in question. 

2.4.1 One particle form factors 

In sharp contrast to the bulk case, in the boundary theory, the boundary operators in 
general may have non trivial one particle form factors (IPFF). Since the multi-particle 
form factors are recursively determined, the IPFF-s are very important inputs to these 
recursions, and their determination is necessarily the first step. The equations for the 
IPFF read: 

F^{d) = R{e)F,{-e) ■ F,{t7r + e)=R{-e)F^{i7T~d), (4) 

where the refiection amplitude R{9) is analytic in the physical strip < Qm{9) < ixjl 
(apart from the presence of finitely many discrete poles on the imaginary axis), and from 
general considerations using the reduction formulae we know that Fx(Q') is analytic for 
< 3m(6') < vr. Note that if Fx{B) is a solution of (jH) then Fx{Q)^{Q) is also a solution 
provided 

^(0) = ^(-0), ^{m ^e) = ^{%Tx -e), 

i.e. if ^ is even and 1m periodic. Therefore one can take "^{Q) = 'ipijj) with y = + e^^ . 

To construct solutions to (jH) we reduce them to a problem already solved in the bulk 
form factor bootstrap. To this end we write -^1(6*) = gi{6)g2{i'n — 6) and suppose that 

g^{e) = R{e)g^{-e) ■ gi{z7i + 9)=g,{zn-e), (5) 

which are nothing else but the bulk two particle form factor equations [H], where the 
reflection amplitude, R{0), plays the role of the S-matrix. Furthermore, plugging this 
product form Fi into (jH) reveals, that g2 must also solve Thus a solution to (jH) can 
be constructed as 

F^i9)=gie)g{^n-9), 

where g{6) is an appropriate solution of (jSj). 

To obtain a solution of ^ we use the following theorem [1^]. If the function h{6) 
is meromorphic in the physical strip < 3m(6') < vr with possible poles at iai, . . . ,iai 
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and zeros at iPi, . . . ,iPk and grows as at most a polynomial in exp(|^|) for \^e9\ oo, 
furthermore it satisfies 

h{9) = M{9)h{-9); M(^)=exp|^ sinh (^^^ | ; h{i7T - 9) = h{in + 9); 

then it is uniquely defined up to normalization as 

_ uu ^i^h iuo - ^inh iuo + ^m [ r dtntf^^^^Sl 

ni=i sinh sinh + [Jo sinhx 

Since the reflection amplitudes are usually expressed as products of the blocks (xj) 

R{0) = n(x.); (x.) = "^^'Hy (6) 

smh(| — i^) 

to use this theorem we need the integral representation of one single block (x),: 

r /"^ rftsinhtfl -x) , {t9\ 

-ix) = exp <^ 2 / ^- sinh — 

^ ^ I Jo t smht / 

Then, if R{9) consists of an even number of blocks, the minimal solution (with no zeroes 
and poles) to eq. (jH) can be written as 



If R{9) contains an extra minus sign, or is the product of an odd number of blocks, 
R = — n(~(^j))i then the factor g{9) necessarily contains a zero at the origin which 
is implemented by putting an extra sinh | into it; thus in this case r^[n{9) = sinh6'r^ijj(6'). 
In the following an important role is played by the appropriate modification of the minimal 
IPFF denoted by r{9) 

r{9) = r^in{9) x — — 
poles 

where the last factor denotes an appropriate number of zeroes and poles at the right places 
(usually the same as in R{9)). 

Thus the general solution of eq.(I3I) can be written as 

F,{9)=r{9)Q,{y), y = e' + e-', 

where the choice of Qi{y) is restricted by the analyticity and the possible asymptotics of 
Fi. It is the Qiiy) in the IPFF that carries the dependence on the boundary operator O. 
Note in particular that if Qi{y) corresponds to the operator O then Qi{y) ~ y^Qiiv) with 
integer N > 1, describes the IPFF of the operator O. 
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2.4.2 Two-particle form factors 

The next step is to investigate the two-particle form factors (2PFF). The novel feature 
compared to the IPFF is that their equations contain also the bulk S-matrix. It is worth- 
while to go through the analysis in some detail since it is straightforward to write down 
the general form of the n-particle form factors once that of the 2PFF-S is obtained. The 
equations for the 2PFF-S have the form 

i^2(^l,^2) = m-^2)i^2(^2,^l), (a) F2(^i,e2)=i?(^2)i^2(^l,-^2) (&) 

F2itn + 61,62) = R{-6{)F2{m- 61,62). (c) (7) 

Note that if ^2(6*1, 62) is a solution to these equations then so is -^2(6*1, 62)H{6i, 62) provided 
if is a symmetric, even and 2m periodic function. 
To construct solutions to eq.Q we write 

i^2(^l,^2) = /(^^l-^2)^(^^l,^2) 

where f{6) is the minimal bulk two particle form factor |14| . i.e. the minimal solution of 

f{6) = S{6)f{-6), f{^,r: + 6) = f{m~6). 

Plugging this F2 into (jlK) reveals that ^ must be symmetric \1/ (6*1, 6*2) = ^{62,61). The 
most convenient way to satisfy ^}p) is that \l/ has the form 

O2) = f{6i + 62)r{6,)r{02)H0i, 02) 

where $ is symmetric and even ^(61, 62) = $(6*1, -62). Finally this F2 satisfies eq.jZt) also 
if $(27r — 61, 62) = $(?7r + 61, 62). The conditions on $ can be satisfied simply by writing 
^(6*1, 62) = 0(1/1, 1/2) where is a symmetric function of the i/i-s {yi = e^* + e~^% i = 1, 2). 
Thus the general form of the 2PFF, compatible with eq.Q is 

i^2(^i, ^2) =r(^i)r(^2)/(^^i-^^2)/(^i + ^2)0(1/1, 1/2), 0(1/1,1/2) =0(1/2,1/1). 

Different choices of the boundary operator O correspond to diff'erent functions 0(i/i, 1/2) in 
this expression. 

2.4.3 Multi-particle form factors 

From the explicit form of the 2PFF it is clear that the general form of the multi-particle 
form factors can be written in the following form: 

n 

F^{6u 62,..., 6^) = Gn{6i, 62,..., 6n) n ri0^) H ^^^^ " ^^^^^^^^ + ^^^^ (8) 

i=l i<j 

where f{6) is the minimal bulk two particle form factor. As a consequence of the form 
factor equations Gn is a 27ri periodic, symmetric and even function of the rapidities: 61, 
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i.e. it is symmetric in the variable yi = 2cosh6'j. When the bulk S-matrix is nontrivial, 
the bulk kinematical singularity equations 



-t Res F„+2(^' + ^^T, 6,61,..., On) = (1 - TT 5(0 - e,)S{9 + %))F„(^^i, . . . , 0„) (9) 



i=l 



give recursive relations linking Gn to Gn+2- (Note that these singularities are absent in the 
two particle case). The advantage of using the i/i-s becomes clear if one tries to describe 
the bulk kinematical singularities: since y{i7r + 6) = —y{0), thus including a (symmetric) 
factor yi + yj in the denominator automatically accounts for the pole. Therefore in the 
following we put 

fa a a \ _ Q™(l/i, I/2 ■ ■ ■ .Vn) 

i<j 

(with Qn being a symmetric function of yi, . . . , ?/„) and then eq.ijHl) give recursive relations 
between the functions Qn- Clearly the actual form of these recursive relations varies from 
model to model since they depend on the bulk S-matrix. The form of the recursions 
depends also on the choice of the IPFF r{6); it is useful to choose an r{6) which gives the 
simplest possible recursion. Writing the 2PFF in the same form as the n-particle one 

^ Q2iyi,y2) 

(p{yi,y2) = ■ 

yi + y2 

then the absence of kinematical singularities requires Q2{y, —y) = 0. 

If the bulk S-matrix is nontrivial and the reflection factor has a pole at y then the 
form factors with odd and even particle number are connected by the boundary kinematical 
singularity equation: 

n 

-I Res F^^,[Q + f , 01, . . . , 0„) = I (1 - n ^(f - • • • > ^n) 

1=1 

The corresponding pole in the n particle form factor can be included as 

Qn{y\,y2-- -.yn) 



Gnidl, O2, ■ ■ ■ , On) 



Uiy^ Uiyi + yj)' 



and the boundary kinematical singularity equation relates Qn to Qn+i- 

The even and the odd particle form factors are also related if the bulk S-matrix has a 
"self fusing" pole describing the 2 particle 1 particle process, which parallels the bulk 
situation (this happens e.g. in the Lee- Yang model). (In this case it is customary to include 
this pole also in f{0)). Since the fusing angle in this process is necessarily 27r/3, one finds 
from bootstrap that in this case the dynamical singularities imply 

-I Res Fn+2{e' + '-^,9-'-^,9i,...,9n)= TFn+l{9, 9i,..., 9n), (10) 
e =e o 6 
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where T is related to the residue of the S-matrix at the self fusing pole: —uesQ^2n±S{9) = F^. 

An important restriction follows on the form factor functions from requiring a power 
law bounded ultraviolet behaviour for the two point correlator of the boundary operators 
(0|O(r)O(0)|0): the growth of the function Fn{6i, . . . ,6n) must be bounded by some ex- 
ponential of the rapidity as ^ — *• oo (i.e. the form factors only grow polynomially with 
particle energy). This can be shown using an argument identical to that in the bulk case 
[T6] . If r (6*) and / (9) are specified in a way to include all poles induced by the dynamics 
of the model, then it follows that the functions Q„ must be polynomials of the There- 
fore in the following we only look for explicit polynomial solutions of the various recursion 
equations. This is a posteriori confirmed since we find as many polynomial solution of the 
boundary form factor equation as many independent local operator exist in the theories. 

2.5 Two-point function 

Once an appropriate solution of the form factor axioms is found it can be used to describe 
correlators of boundary operators. The two-point function of the boundary operator O 
can be computed by inserting a complete set of states 

oo 

(0|C»(t)C»(0)|0) = J2 TTT^ / d9,d92 . . .d^„e-^'"*S»^°^^'''F„F+ (11) 

n=o ^ Jei>e2>- >e„>o 

where time translation invariance was used and the form factors are 
F„ = (010(0)1^1, 92,..., 9^),^ = F^{9u 92,..., 9^) 

and 

= in{9l, 92,..., 9n\OmO) = (ZTT + ZTT + 9^_,, ...,171 + 9,) 

which, for unitary theories, is the complex conjugate of the previous one: = F*. In the 
Euclidean (r = it) version of the theories the form factor expansion of the correlator for 
large separations converges rapidly since multi-particle terms are exponentially suppressed. 

The identification between solutions of the form factor equations and operators of the 
theory is a central issue. One possible way is to analyze the behaviour of the boundary 
correlators for short distances. Although on general grounds one may expect the form factor 
expansion to converge rapidly only in the infrared (large volume) regime, the examples from 
the various bulk theories, where the form factor expansion converges even in the UV domain 
(see e.g. [l_7j), suggest that similar behaviour may happen in the boundary setting as well. 
If the theory can be described as a relevant perturbation of a conformal field theory, then 
in the UV domain the two-point function must follow a behaviour dictated by this limiting 
theory. The short distance singularity exponent is related to the scaling dimension of the 
operator O and can be calculated from the asymptotic growth of the form factors. 
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3 Model studies 



In this section we carry out a detailed investigation of the solutions of the form factor 
equations in four different models. The first two models (the massive scalar field with 
linear boundary interaction and the Ising model interacting with a boundary magnetic 
field 0) are free in the bulk and the correlation functions are known explicitly, thus the 
form factors obtained from the explicit field theoretic solutions can be compared directly 
to the solutions of the form factor equations. In both cases we find that the space of 
appropriate polynomial solutions of the FF equations can be identified with the space of 
local boundary operators obtained from the explicit construction. In the case of the Ising 
model we also show how the conformal dimensions of the various operators of the ultraviolet 
limiting BCFT can be obtained from the solutions of the FF equations. The third and 
fourth models, namely the scaling Lee- Yang and the sinh-Gordon models with integrability 
preserving boundaries are among the simplest boundary integrable theories. In contrast 
to the previous cases they cannot be solved directly so one has to rely upon the solution 
of the form factor equations. Since these models contain nontrivial bulk interactions the 
recursion relations connecting the multi-particle form factors are no longer trivial, and in 
these cases we investigate their solutions in detail. 

3.1 Massive scalar with linear boundary interaction 
3.1.1 Direct calculation 

The free massive scalar field t) restricted to the negative half-line x < subject to 
linear boundary condition 

d:,^x,t)U=o = -HHO,t)-%). (12) 
can be described by the following Lagrangian: 

/I 1 m"^ \ \ 

This one parameter family of linear boundary condition interpolates between Neumann 
dx^\x=o = (for A = 0) and Dirichlet $|x=o = "^"0 (for A — > 00) boundary conditions. 
Since for any A we are dealing with a free theory it can be solved explicitly. The mode 
decomposition of the field is 

$(x,t) = Ae'""+ [ 47T|«(^)e~'"^''^*(e'''"+^(fc)e"'''")+a+(A;)e*"('=)*(e-''="+i?(-A;)e*'=")| 
Jo uj{k) I J 

where A = ^-r*&o and 

m+A 

m = ^ 

^ ^ k + tX 
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is the reflection factor on the boundary at x = 0. The creation/annihilation operators are 
normalized as 

[aik),a+{k')] = 27iujik)S{k~k') , k , k' > 
The boundary two-point function can be calculated easily 

(0|<l>(0, t)<l>(0, t') |0) =A^+ _^^^-Mm-t ) (1 ^ (1 + 

Jo Z-nuj^k) 

By comparing this expression to the form factor expansion of the two-point function (fTT| . 
the form factor of the elementary field can be extracted: 

(0|$(0,t)|^) = e-'"('=)* (1 + i?(A;)) 

The same result can be obtained by taking the general (space-dependent) two point function 

(0|r(^<l>(x,t)$(x',t')) |0) = A2e'"("+"') 
and using the boundary reduction formula IH] 

/O /* oo 

rfx' / dt e'"(^)*' cos(p(^)x') js^, - dl + + | 
oo J — oo 



(0|T(^$(x,t)$(x,t)j |0) 

where uj{9) = mcosh(^) and p(6') = msinh(^). Performing explicitly the calculation yields 

(0|<l>(x,t)|^) = e-^^(^)*(e'P(^)" + i?(e)e-**'(^)") (13) 
which for the form factor of the operator $(x = 0, t) reads as 

{omo,t)\e) = e-'^^'^' {1 + R{e)) 

Introducing r = —it one also gets 

{o\d!^<i>{o, o)\e) = uj{eY (i + R{e)) , n>o. 

Clearly these operators have no multi-particle matrix elements. It is important to realize 
that dx^{0,0) is not an independent operator since the boundary condition eq. (fT2ll relates 
it to $(0,0), thus the set of independent boundary operators having only one particle 
matrix elements is given by (9"$(0,0). To obtain multi-particle matrix elements one has 
to consider 

(0| : $(xi,ti)...<l>(xfc,4) : \9^...ek). 
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Using the analogous boundary reduction formula and the Wick theorem we obtain 

(0| : : 1^1... ^fc) = (14) 

where the ellipses at the end represent additional terms which make it completely sym- 
metric in all coordinates. From this expression one can extract the form factor of the most 
general boundary operator of the theory 

(0| : d'^^'^O, 0) . . . dl'<!>{0, 0)■.\9^... 9^) = u;{9,r (1 + ^(^i)) • • • ^iO^r (1 + ^(^fe)) + • ■ ■ 

where again a complete symmetrization in the 9i rapidities is understood. Checking the 
leading asymptotic behaviour of these form factors gives that for all 9i ^ 9 large they 
grow as e^^, where N = ni + ■ — h is the total number of derivatives in the expression. 
We note that we have as many operators for a given as many partition N has into the 
numbers 1, 2, . . . , /c. This can be seen by writing N = Ni + 2N2 + ■ ■ ■ + kNk and associating 
to it the operator with rii = Ni + N2 ■ ■ ■ + Nk, n2 = N2 + ■ ■ ■ + Nk . . . Uk = Nk derivatives. 

The Dirichlet boundary condition {R = —1) can be obtained in the A ^ 00 limit. 
Clearly $|x=o = $0 is a c- number and the Dirichlet boundary condition does not connect 
the operator dx^\x=o to $|a;=o- We can extract, however, the form factors of the operator 
9a;$(0, t) from that of $(0, t) by taking the A ^ cx) limit carefully in ^H^ : 

{0\d^<l>{0,t)\9) = e-'^^'^Hip{9) 

and for its derivatives 

(0|(9;^9^.<l>(0,0)|^) = u{9Y2ip{9), n > 0. 

This can be extended similarly to the most general operator as 

(0| : d^.'dM^, 0) . . . 9,"'=9,.$(0, Q):\9r...9k)= u{9,r^2zp{9,) . . . uj{9kr'^2ip{9k) + • • • 

where again a complete symmetrization in the 9i rapidities is understood. Checking the 
leading asymptotic behaviour of these form factors gives that for all 9i 9 large they grow 
as e^^, where N = k + ni + ■■■ + Uk is the total number of derivatives in the expression. 

3.1.2 Solving the form factor equations 

The bulk S-matrix of the theory together with the reflection factor are 

, , sinh^-i^ f B\ f B\ A nB 

S(9) = 1, R(9) = ^ = - 1 + _ , - = sin , 

smh9 + i^ V 2 J \ 2 J ' m 2' 

where the block notation (jH)) is used to express R{9). As a consequence of 5* = 1 the 
minimal bulk form factor is trivial: f{9) = 1. To determine the IPFF we note that this 
reflection factor is identical to the two particle S-matrix of the sinh-Gordon model if the 
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above identification of parameters is done. Therefore eq.(ini) in this case is identical to the 
equation for the minimal bulk form factor of the sinh-Gordon model. Choosing for g{6) 
the solution given in jH] (described in detail in the sinh-Gordon section) gives 

r{e) = 2g{e)g{zn -6) = -^^^^ = 1 + R{e) 

smh (J + 1— 

rn 

Clearly this corresponds to the form factor of the operator $(0, 0) — A. 

Now we demonstrate that the number of independent solutions of the form factor 
equations coincides with the number of local boundary operators. In this case the general 
Ansatz (jH} takes the following form 



On) = PniVl, ■ ■ ■,yn)Ylr{ei 



where i/i = 2 cosh 9i as before. Since the bulk S'-matrix is trivial there are no bulk/boundary 
kinematical singularities and P„ is a completely symmetric polynomial in the i/i-s. One can 
count the independent solutions of the BFF equations by counting the possible solutions 
for P„. If Pn has degree N then the solutions are given by the partitions of N into the 
numbers 1, 2, ... n in the following way. Since the completely symmetric polynomials of n 
variable are generated by the (Tj-s (elementary symmetric polynomials of degree i) one can 
write: 

n n 
i=l i=l 

It is clear that this space has the same dimension as the space of boundary operators 
having only n particle matrix elements with asymptotic growth e^^ . 

The Dirichlet (A oo) limit for the operator 9j,$(0, 0) can be obtained using its 
relation to $(0, 0) via the boundary condition (fT^ as we did in the Lagrangian framework. 

3.2 Ising model with boundary magnetic field 
3.2.1 Direct calculation 

The Ising model with a boundary magnetic field can be described by a free massive Majo- 
rana fermion perturbed at the boundary 0. In Minkowskian formalism the Dirac equation 
can be obtained form the Lagrangian: 

The gamma matrices are chosen as 
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in order for the Dirac equation to be real: 

- dt -m J 

Thus the Majorana condition corresponds to taking ^' real. Using the component notation 
the Lagrangian of the boundary field theory takes the form 



C = -iQ(^-x) (^^^+{dt - d:c)'ip+ - ^^-{dt + d^)ip- - mil^+ip-^ - id{x)UB 

where ^ ^ ^ 

Ub = -V'+V'- + ^aa + -ha('ijj+ + -0-) 

The operator a is a boundary fermion = 1, which implements the boundary condition 

Since the theory is free we can solve it explicitly. The mode expansion of the fermionic 
fields are 

where u±{9) = v^{9) = ^Jme^^^ are the spinor amplitudes, R{&) is nothing else but the 
one-particle refiection factor 

. -, / ^TT Q\ sinh^ + m ^ 'b? 

R(e) = i tanh — — , k = 1 

^ ^ V 4 2 J sinh O-m ' 2m 

and the creation/annihilation operators are normalized as 

{b{e),b+{e')} = 27rS{e-e') 

The boundary two point function can be calculated explicitly: 

(o|v^±(o,t)v^+(o,t')|o) = J -e— ) {u^{e) + R{e)u±{-e)) k(^) + R{-e)v^{-e)) 

(15) 
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which, as compared to the form factor expansion gives 

(0|^±(0,t)|^) = e-^-(^)*(«±W + ^(%±(-^)) (16) 

These two operators are not independent since they are related by the boundary condition 
and so there is only one boundary fermion field, say V+- As a result, the algebraically 
independent operators at the boundary are the fermion field and its derivatives d!^ip+- 
Note that drcip+\x=o is not an independent field, as it is determined by the Dirac equation 
in terms of drip+\x=o and iIj_\x=o. As a consequence of the fermionic nature of the field 
the most general boundary operator has the form d!^'^ip+d'!^^ip+ . . . d'!^'"ip+ where ni > n2 > 
■■■ > Uk (the inequalities are strict, in contrast to the bosonic case discussed earlier). 
N = rii + n2 + ■ ■ ■ + Uk is called the level of the operator, and operators at level can 
be brought in one-to-one correspondence with partitions of N into the numbers 1, 2, . . . , /c. 
For a partition 

N = kNk + {k- l)Nk-i + --- + 2N2 + Ni 

we associate the operator above with Uk = Nk, Uk-i = Nk + Nk-i] . . . ; rii = Nk + Nk-i + 
■■■ + N^. 

(Had we included also the operator a we would have had to perform a GSO type pro- 
jection, leaving only non-fermionic operators. This would amount to keeping all operators 
with an even number of fermion factors plus all odd ones multiplied with a factor a, but 
this would lead to the same number of operators.) 

3.2.2 Solution of the FF bootstrap 

Using again the block notation ^ the S-matrix of the theory and the reflection factor are 

[5] 

3(6) = -1 , R,(e) = [x] (^-^^ , [x] = ix)il-x) 
where x is related to the boundary magnetic field as 

sinvrx = 1 = k 

2m 

For h = we recover the free boundary condition with 

-Rfree(^) = 

which has a pole at ?| corresponding to the fact that the ground state is doubly degenerate. 
In contrast to the generic situation the pole at i| is a dynamical pole and not a kinematical 
one (since the field has no vacuum expectation value at all). The h ^ oo limit corresponds 
to the fixed boundary condition (when the Ising spin takes a fixed value at the boundary), 
and the refiection factor is 

-Rfixed(^) = 
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which has no pole in the physical strip at all. The minimal one particle form factor for the 
fixed case can be calculated directly using the recipe in Section 2 



'^'^^'^^'^^ sinh(f + ^f) 
For the free case we include the dynamical singularity into the IPFF 

^.sinh^ . , 

The simplest solution which interpolates between these two cases and has a pole exactly 
at the location of the boundary dynamical singularity of the reflection factor is 

, . sinh^ cosh6' + z(l — k) 

''^ ' " sinh(f + zf) sinh^-z/€ ' ^ 

This expression is the same we obtained from the exact solution of the model (fTKll . The 
minimal bulk two particle form factor is simply 

/(^)=sinh^. 

Since R{0) has no kinematical pole at boundary kinematical singularities are absent, 
and since the bulk S'-matrix is —1 there are no bulk kinematical singularities either. Thus 
we look for the n particle form factors in the form 

Fni9i, ...,9n)= Pnivi, . . . , y„) n ^(^') n - ^^•)/(^* + ^^•)' 

i i<j 

where i/i = 2 cosh and P„ is a completely symmetric polynomial in the i/i-s. Taking into 
account the special form of f{9) the form factor simplifies to 

. . . , = Pnivi, . . . , n ^(^^) Hiy^ - yj) (18) 

i i<j 

The independent solutions are counted in the same way as in the bosonic case, i.e. by the 
partitions of into the numbers 1, 2, . . . , n and are generated by the cxj-s. It is clear that 
the dimension of the space they span is the same as the one of the boundary operators 
obtained from the direct calculation. 

Since the UV limit of this theory is a boundary conformal field theory one can go 
further than in the bosonic case and calculate the UV dimension of the various boundary 
operators. As the form factor equations are not coupled we can choose a basis among 
operators consisting of those having matrix elements only with a certain fixed number of 
particles. 
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If the operator has only one-particle matrix element then its correlator is 

(0|C(r)O(0)|0) = / i^|FC?(^)|V"^^°^^'" (19) 
Jo 27r 

where F^{0) = r{6)Pi{y). Plugging (|T7|l into (fm we obtain the exact correlator (fTKll . If 
the operator O goes to a scaling operator in the UV limit (r — > 0) then the exact correlator 
has the short distance asymptotics r~^^, where A is the appropriate scaling dimension in 
the ultraviolet BCFT. In (fT9| the singularity comes from the large 9 asymptotics of the 
form factor. If |Ff (6')p diverges as here, then the corresponding weight is A = |. 
Taking the simplest solution Pi{y) = 1 the weight is A = | which corresponds to the 
boundary fermion field. Choosing Pi{y) = cr^iy) corresponds to the n-th derivative of this 
operator which has weight ^ + |. 

Similarly we can analyze an operator having n-particle matrix element only. The cor- 
responding correlator is 

(0|OH(r)O(.)(0)|0) = f^---f ^Ijl^niOu . . . , V-(-^^^^+-+-^^^")- 

The operator which has the mildest UV behaviour corresponds to P„ = 1. The correspond- 
ing form factor square for large 9-s behaves as 

\F^i9,, . . . , cx exp{9in + n{n - 1)) = e'^\ 

2 

thus the UV dimension of 0(n) is A = The explicit boundary operator which has 
nonzero matrix elements only with n particle states and has the mildest UV behaviour is 

with dimension A = | + ; therefore it can be associated to 

To match the descendent operators, note that to any partition of = ki+2k2+- ■ ■+n-kn 
there exists a solution of the form factor equations with Pj^ = . . . cr^" . The number 
of such polynomials is the same as the number of descendants of 0(„) at level A^: to the 
given partition we can associate the operator 

Conversely, given an operator of the form 



„2 



of weight + one can define a partition as kn = Pi, ^n-i = P2 — — 1, • • • and thus 
associate a polynomial solution of the form factor equations with appropriate asymptotic 
behaviour. It is important to emphasize that we do not claim that the form factor related 
to P^ belongs to the operator above, what we have shown is only that the dimension of 
the space of operators with certain scaling dimension is the same as the dimension of the 
solution of the form factor equations. 
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3.3 The boundary scaling Lee- Yang model 

The scaling Lee- Yang model with boundary is a combined bulk and boundary perturbation 
of the boundary version of the M.{2/b) minimal model [HIIHl- In the bulk the perturbation 
is given by the unique relevant spinless field 0, at the boundary the perturbation depends 
on which of the two possible conformal boundary conditions is present in the unperturbed 
model. One, denoted by 1 in [TH], does not have any relevant boundary fields - thus can 
have no boundary perturbation either -, while the other, denoted $ in [TH], has a single 
relevant boundary field ip with scaling dimension —1/5. In this latter case the general 
perturbed action is 



A 



.4# + A dy dx(f){x,y) + h / dy(p{y) 



where Aq, denotes the action for A4{2/5) with the $ boundary condition imposed at x = 0, 
and A and h denote the bulk and boundary couplings respectively. The action of A\^i is 
similar, but the last term on the right hand side is missing. If A > then in all cases the 
bulk behaviour is described by an integrable massive theory having only a single particle 
type with the following S matrix 



S{9) 



X, 



_ sinh (I + ^) 
- sinh (I - ^) ■ 



The pole at ^ = ^ corresponds to the "ip^ property", i.e. the particle appears as a bound 
state of itself. The minimal bulk two particle form factor which has only a zero at 6' = 
and a pole at 6* = ^ in the strip < Qm{9) < n has the form [T7j : 



m 



y 



y + l 



-v[tn 



-in 



y = e" + e 



where 



v{d) = exp <^ 2 



T 



. et sinh I sinh | sinh | 

e IT — ^ " 



sinh t 



In the boundary theory with the perturbed $ boundary, the reflection amplitude of the 
particle depends on the strength of the coupling constant of the boundary held as E 



R{9)^ = Ro{9)R{b, 



b + 1 



6-1 



where the dimensionless parameter b is related to the dimensionful h as 



h{b) = sin ( {b+^)^)miXf/^h 



criti 



''crit 



3 4 1 

_7r52554 



sm 



2lT 



r 



r(i) 
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and m(A) is the mass of the particle giving the overall scale in the infrared description. In 
the case of the 1 boundary the reflection amplitude is the parameter independent expression 

^<^" ^ (i) (i) (-r 

Note that -R(^)i is identical to -R(6')$ at 6 = and so both have a pole at 6 = i7i/2 coming 
from the (|) block but their g factors differ in a sign [T8] . 

3.3.1 Lee- Yang model with perturbed $ boundary 

We consider flrst the Lee- Yang model with perturbed $ boundary. The IPFF correspond- 
ing to i?(6')$ is chosen as 

i sinh e 

r[0) = T ^ u{9), 

(sinh 6 — ism ^ ' ) (sinh 6* — i sin ^ ' ) 



where 



u{9) = exp 



dt 
T 



2 cosh - cos 



sinh f 2 



2 / vr 



sinh 7^ + sinh | — sinh | 



sinh^ t 



Note that r ~ 1 at y — > cxd, and r(6') satisfies the r{6 + zvr) = r(6')* reality condition for 
real 9. The general n-particle form factors have the form 

, I j + * 

where we separated a normalizing factor if„ from the polynomials Qn- The various F„-s 
are related to each other by both the kinematical and the dynamical singularity equations, 
since the S-matrix is nontrivial and also has a ip^ pole with V = z22 34. In addition, these 
Fn-s also have to satisfy the equation coming from the residue of the pole at 6' = i7r/2: 

-2 Res F„+i(^, 9,,...,9r,) = ^[l-f[ ^(^ - 0,) ] F^, . . . , (21) 

\ i=i / 

Our strategy is to solve the recursion equations coming from the first two conditions first 
and check whether the solutions also satisfy the third requirement (f2T| . By choosing the 
normalizing factors and introducing the useful quantities (3k 

^— (3k{h) = 2cos-{h + k), keZ, (22) 

22t>(0) J t) 

the overall normalization N drops out and the recursion equations coming from the dy- 
namical (resp. kinematical) singularities read 

Q2{y+,y-) = {y'-P'^s)Qi{y), 

n 

Qn+2{y+,y-,yi,- ■ - ^Vn) = Qn+i{y,yu---,yn){y^ - (3l3)Y[{y + yi), n > 0; (23) 

i=l 
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where 



Qn+2{-y,y,yi,---,yn) = Qniyi, ■ ■ ■ ,yn) {y^ - (3^i){y^ - Pl) Pn, 



Yliyi - y-){yi + y+) - Wiyi + y-){yi - y+) 



2(^/4 



and 



y^ = ujx + uj ^] 



y-] 



y- 



1=1 



i=l 



X 



6 3, y = X + X 



-1 



(24) 
(25) 

(26) 



Next we present the minimal solution of these recursion equations up to n = 3. The 
solution is called minimal, if the leading overall degrees of the F„-s in all of the y variables 
are the smallest possible ones. Of course the solution also depends on the input function 
Qi(j/i). Since Fi can have no pole at 6* = i7r/2 while r{6)/y has one, Qi must be chosen to 
cancel this pole; the choice with the minimal degree is Qi{yi) = yi = ai{yi). Using this as 
input, we find from (OTH 1^ the unique solution 



Qiiyi 



0-1, 



3;) 



Qsiyi, 2/2, 2/3) = en [cri(cr2 + P-i){cr2 + /5i) - {^2 + S){aia2 - 0-3 



(27) 



The remarkable property of this solution is that it contains no free parameters. A simple 
counting of the various powers shows that the leading overall degree of Fi, F2, and F3 
vanish. 

To check eq. (f2T| we need the following relations following from the explicit solution 
(f27|) and from the various identities among the /3fc-s: 



^2(0,2/2) 

Q3(0,?/2,l/3) 



o-i(y2)(3-/3! 



3J^ 



(3i(3-i(Ti{y2,y3)Q2{y2,y3) = - P-3)My2,y3)Q2{y2,y3)- 

Indeed using them in eq. (OT|l leads to a consistency condition on the ratio of the Hn-s: 

Hn+l 



^ -^(9(3-/3^3) 



n 



1,2. 



Since 



ru- 



4M(if) 



and g^ = t2i3y/\2-V3)'/' ^^^-' 



z ^(v^-/5_3)2' 73-/3-3' 

the b dependence cancels from the consistency condition and using the actual form of the 
Hn-s leads to 

^i^ = v^(2-v^)V2, 
which we checked numerically up to 7 digits. 



24 



To test these form factors numerically against the predictions of conformal field theory, 
we take the spectral representation of the boundary two-point function 

°° poo JQ JQ 

(o|o(t)0(o)|o) = Y.j^ 

which we truncate to the first few terms in the boundary form factor expansion. Since the 
minimal solution of the form factor problem has the mildest UV behaviour it is natural to 
assume, that in the UV it corresponds to the boundary field (p. Therefore {0\O(t)O{0)\0) 
obtained from the FF expansion must be compared to the short distance expansion: 

{0\m^-(p{t)m^(p{0)\0) = -(mt)i + (mt)ic^^(m^(^) + . . . 

where appropriate powers of m were inserted to make the expression dimensionless and 
the fusion coefficient is 

— 



'n 


'1- 


)r( 


;§) 




'3' 
^5. 


)r( 





2 

while the (6-dependent, dimensionless) vacuum expectation value 



5 cos(f) 
6K„t cos(f^(26+l)) 



is given explicitly in [TH]. In analogy with the bulk case |7| we choose the normalization 
factor N in (f^ as the vacuum expectation value of the boundary field 

N = (m^if) (28) 

With this choice the boundary form factor expansion gives 

{0\O{t)O{0)\0) = \F^\^- / ^\F^\\-^i-oshe 

Jo 27r 



d6id62 I 0^ /I \|2^-mt(coshei+coshe2) 



6(27r)3 

The two expansions are compared on the next figure 
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0.007 0.008 



0.009 



0.01 



where b = —1.05 and the dimensionless correlation function is plotted against mt. The 
predicted UV behaviour is given by the continuous line and the numerically determined 
form factor expansion truncated at 1, 2 and 3 particle intermediate states is denoted by 
the symbols □, x and o, respectively. 

We checked that the agreement between the form factor expansion truncated at three 
particles and the UV CFT prediction holds for various values of the parameter b: indeed 
as we change b the two curves move together. The agreement above also confirms our 
choice (f^ for the normalization of the form factors. Based on all these we associate the 
boundary operator corresponding to the minimal solution of the form factor axioms to the 
one, that in the UV limit becomes the ip field of the boundary Lee- Yang model. 



3.3.2 Lee- Yang model with the 1 boundary 

The IPFF corresponding to the parameter free R{9)i is chosen as 

ri{9) = i sinh 9 u{9), 

where u{9) is the same as in the previous subsection. Note that ri also satisfies ri{9 + i7i) = 
ri{9)* but its asymptotic behaviour ri ~ at ?/ — oo is different from that of the r in the 
previous case. Since R{9)i also has a pole at ^ = i7r/2 we introduce a similar Ansatz as in 



(EHl) 



9n)=KQn{yU---,yn)l[-^l[ 



rii9,)yjf{9,-9,)fi9, + 9,) 



i<j 



Vi + Vj 



(29) 
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with = 4"if„ where Hn is the same as in (f22|l . Then one finds the following recursion 
equations for the Qn-s from the dynamical (resp. kinematical) singularity equations: 

Q2{y+,y-) = Qi{y), 

n 

Qn+2{y+,y-,yi,---,yn) = Qn+i{y,yi, ■ ■ ■ ,yn) Yliy + yi), n>0; (30) 

i=l 

Qn+2i-y, y,yi,..., yn) = Qniyi, ■■■,yn) Pn, (31) 

where the various symbols are the same as in eq. ()23|24l) . Up to n = 4 the unique minimal 
solution of these recursion equations with the input Qi{yi) = cti is 

Qliyi) = (^U <52(2/l,2/2) ~ <53(2/l,2/2,l/3) ~ O-?, <54(yi,2/2,2/3,y4) ~ 0-?(cr2 + 3). 

It is easy to show that the leading overall degree of the first four form factors Fi, . . . , F4 is 
two. This indicates that the operator that corresponds to this set is different from the one 
encountered in the case of the perturbed $ boundary. Therefore in the conformal limit this 
operator is different from the (p field and this is in accord with the fact that only the identity 
operator and its descendents can live on the conformal boundary condition 1. Based on 
the asymptotics of the form factors for large 9 the corresponding operator has ultraviolet 
dimension 2 and can be identified with the unique such operator in the conformal vacuum 
module which is the L_2 descendent of the identity. This identification is further confirmed 
by comparing the numerically obtained truncated form factor expansion to the conformal 
two-point function. 
Since 

rl(^|) = -t^(4), and ^71 = -^2(3)^^(2 - v^)i/2, 

(see also [12]) one can readily show that these four form factors also satisfy the equation 
coming from the residue of the pole at 6^ = i7r/2. 



3.4 The boundary sinh-Gordon model 

The sinh-Gordon theory in the bulk is defined by the Lagrangian^: 



m 



(cosh 6$ - 1) 



It can be considered as the analytic continuation of the sine-Gordon model for imaginary 
coupling (3 = ib. The S-matrix of the model is 



S{0) 



1 



B 
2" 



B 
' 2" 



B 
2" 



B 



2b^ 



87r + 62 



^Note that the parameter b is used here with a different meaning compared to the former case of the 
boundary Lee- Yang model. 
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The minimal bulk two particle form factor belonging to this S-matrix is ^Ij 



dx 



sm 



X 



XilTT 
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sinh ^ sinh(l — gjnh £ 



2 ^ 2 



sinh a; 



and it satisfies 



sinh 6 



sinh ^ + 2 sin ^ 



B 



(32) 



The sinh-Gordon theory can be restricted to the negative half-line, but the integrability is 
maintained only by imposing either the Dirichlet 



D 



$(0,t) = $ 



or the two parameter family of perturbed Neumann 



VemO^t)) = Mo cosh ( -($(0,t) - ^o] 



Mn 



boundary conditions. The latter interpolates between the Neumann and the Dirichlet 
boundary conditions, since for Mq = we recover the Neumann, while for Mq oo the 
Dirichlet boundary condition with = $o- The refiection factor which depends on two 
continuous parameters can be written as 

R{9) = Ro{9)R{E,F,e)=(^ 



B 



B 

-1 

4 



E 



in terms of the parameterization used in |2D]. They are related to the parameters of the 
Lagrangian as 



E F 
cos — (6^ + Svr) cos — (6^ + Stt) 
Id lo 

E F 
sin — (6^ + Stt) sin — (ft^ + Svr) 
lb lb 



Mo , 6$o 
■ cosh 

6$o 



Merit 

Mo 



M 



sinh 



crit 



where Merit = w^rMp~/^- Note that for generic values of the parameters (i? 7^ 0, 

F 7^ 0) this refiection factor has a pole at 6' = i7r/2 coming from the (|) factor. Imposing 
Dirichlet boundary condition instead of the general one corresponds to removing the F 
dependent factor from R{9). Then the remaining parameter E is related to the $0 boundary 
value of the sinh-Gordon field as. E = iSb^o/^b"^ + Svr). 



3.4.1 Sinh-Gordon model with $0 = Dirichlet b.c. 

This case is interesting because E = implies that the pole at ^ = i7r/2 is absent in this 
case. Therefore the equation coming from the residue of this pole is also absent and the 
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form factors are less restricted. The IPFF corresponding to the reflection amphtude on 
the E = Dirichlet boundary is 



where 



u{9, B) = exp 



-2 1 ^ 

X 



ATT X 

cos{—-e)- - 1 

I TT 



cosh I f . , xB . , ,^ B 



X . , X 



sinh h sinh(l ) — h sinh — 

sinh X \ 4 2 2 2 



Note that ro ~ y at ?/ oo and has no pole at 6* = i7r/2. At -B = - which corresponds to 
a free theory - m(6',0) can be integrated explicitly yielding ro(6')|B=o = (— i sinh 6')/2; and 
this, apart from a trivial normalizational phase coincides with the IPFF for a free scalar 
with Dirichlet b.c. (As discussed in Section 3.1.1 in this case dx^iS^.Gi) is the operator 
having one particle matrix element only). 

We write the n particle form factors in the general form: 



yi + 

Since there is no self fusing pole in the S matrix of the sinh-Gordon model, the -F^-s 
are related only by the kinematical singularity equations. Choosing the ratio of the -ff„-s 
appropriately the recursion equations originating from here take the form: 

Qn+2{-y, 1/1, ... , Vn) = -QniVl, ■ ■ ■ , Vn) Pn, 

where P„ is given by eq. ()25|26j) with u = e''"^ . 

As tq has no pole at ^ = z7r/2 , one can have a minimal solution of this recursion 
equations starting with Qi = 1 which has non vanishing form factors for odd particle 
numbers. We calculated up to = 5 and found that the solution is uniquely given by 

Qi{yi) = 1, <53(i/i,z/2,z/3) = -f^i, 

QbiVi, • • • ,2/5) = o"i[o-30-2 - {uj + uj~'^ya5 + {uj- - (a; - cj"^)^(a3 + aida)], 

with all the Fi, F3 and F5 form factors having leading overall degree one. There is a unique 
local operator with this property, namely d^^. 

Of course one can also find non vanishing solutions with even particle numbers also 
starting with a non trivial Q2. Since F2 can have no kinematical singularity, the minimal 
choice is (52(1/1,1/2) = cti- With this input we obtained again a unique solution 

^2(1/1, 1/2) = CTi, Qiiyi, . . . , 1/4) = cr?(c^2 -{uj- ^^^f), 

where both the F2 and the F4 have leading overall degree two. 

For $0 = (-£" = 0) the $ — $ reflection symmetry of the bulk sinh-Gordon model 
survives also in the boundary theory. Therefore the boundary operators can be classified as 
even or odd ones, having only non-vanishing even or odd particle form factors, respectively. 
Thus the second form factor family can be identified with the operator {dx^Y . 
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3.4.2 Sinh-Gordon model with $o 7^ Dirichlet b.c. 

For $0 7^ (i? 7^ 0) the reflection factor acquires a pole at 6 = in/ 2 due to the fact that 
the fleld has a nontrivial vacuum conflguration. At the same time the reflection symmetry 
of the bulk sinh-Gordon model is violated in the boundary theory. Therefore the boundary 
operators cannot be classified into representations of this symmetry, and the equation 
coming from the residue of the pole at 6 = i7r/2 connects the form factors with even and 
odd particle numbers. Note that now this equation plays an essential role as it is the only 
one that relates these two sets of form factors to each other. 

The IPFF corresponding to the refiection amplitude on the E ^ Dirichlet boundary 

is 

smh — 2 sm 7 2 

where u{9, B) is the same as in the previous case. Note that r^: ~ y at y ^ 00. 
Writing the n particle form factors in the general form 

and choosing the ratio of the Hn-s appropriately the recursion equations originating from 
the kinematical singularity equation take the form: 

Qn+2(-l/,l/,l/l, . . . ,1/n) = (2/^ - 4C0S^7)Q„(?/1, . . • , 1/n) Pn, (33) 

where P„ is given by eq. ()25f26|) with uj = e™^ . Next we show how the equation coming 
from the residue of the pole at 6* = ot/2: 

-I Res Fn+i{e, ^1, . . . , = ^ I 1 - n 5(4 - d,) ) F„(0i, ....On), (34) 

helps to eliminate the arbitrariness in the minimal solution of the recursion equations. 

Qiijji) = 0"! is the minimal choice that guarantees that Fi has no pole at = i7r/2. 
Using this in the recursion equation (j33|) gives that the most general has the form: 

^3(1/1,1/2,2/3) = -Cr?(0-2 + 4C0S^7) + {A + 5(7i)(cri(72 " (^3), 

where A and B are arbitrary constants. Eq. (|r?Hl leads to the following relation between Q3 
and Q2: 

H3rE{'i-)Q3{0, y2, 2/3) = 9e2 sin ^i^2(Xi(?/2, 2/3)Q2(2/2, 2/3)- 

Since 

(53(0,1/2,2/3) = {-0-1(0-2 + 4cos^7) + {A + 50-1)0-2} , 

the expression in the curly bracket should be proportional to Q2- This observation fixes the 
values of A and B: Q2 has to be proportional to ai to guarantee that F2 has no kinematical 
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singularity and this requirement is met only if A = 0, while F2 has a leading degree not 
exceeding that of Fi and F3 if 5 = 1. Thus with these two requirements one obtains a 
parameter free solution starting with Qi = up to n = 4 it has the form: 

(52(2/1, 2/2) 4 cos^ 70-1, Qsivi, 2/2, ys) = -0-10-3 - 4 cos^ 7cri, 

nB 

QiiVi, • • • , 2/4) ~ -4 cos^ 7(0-1 era + 4 cos^ 70"^) (as + 4 sin^ — ). 

Note that both Q2 and Q4 vanish for 7 = — 7r/2 (E = 0). Furthermore for £" = one also 
has Qi/ui = 1, and Qzl {yiy2y-i) = — o"!, thus the solution goes over smoothly into the one 
with E = Since 

.TT, 1 , 2(1 + cos ^ + sin ^) COS7 



''"E{i-) = : , and gE 



2 1 — sin 7 ttB 1 — sin 7 

V ~ 

the 7 dependence cancels from the ratios of H^/Hi and H4/H2 when we use eq. (jr?ill for 
n = 1,2,3: 

Hs f ttB nB\^ TvB 

—i-r- = —i-rr = 1 + cos — — h sin — — 4 sin ■ 



H2 Hi \ 4 4 y 2 

In the 6 — i> limit these ratios vanish, therefore the higher form factors decouple in accord 
with the fact that the kinematical singularity axiom becomes trivial for the free field theory. 



4 Conclusion 

In this paper we treated the form factor bootstrap for boundary operators in integrable 
boundary quantum field theory. Although there have been earlier treatment of form factors 
for specific (mainly lattice) models [IHl IHl 1121 , none of these has actually given a complete 
formulation similar to the axiomatic approach by Smirnov for the bulk case f^. The present 
work initiates an extension of this axiomatic program to boundary fields. 

We have given a complete axiomatization of the properties of boundary form factors, 
derived from first principles of quantum field theory (unitarity and the boundary extension 
of the LSZ reduction formulae). In particular, the axiom describing boundary kinematic 
singularities is an entirely new result of this paper, as this has never been treated before 
in any previous study. We have shown that these axioms are consistent with many known 
aspects of integrable boundary field theory. In particular, the relation between the residue 
of the reflection factor at i'7i/2 and the one-particle contribution to the boundary state, 
noted previously in the context of finite size effects, was confirmed once more as a necessary 
condition for the consistency between the boundary and the bulk kinematical axiom (the 
only exception to this relation is when the bulk is free, but then the two axioms are trivial). 
Therefore it seems that this particular relation is a consequence of integrability and the 
existence of a nontrivial bulk scattering matrix. 

We then proceeded to give a systematic method to solve the boundary form factor 
axioms for the case of diagonal scattering. The solution is a natural generalization of the 
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bulk case, but necessitates the introduction of a minimal boundary form factor function in 
addition to the already known minimal bulk form factor. The periodicity, permutation and 
reflection axioms can then be solved by a general Ansatz, and the residue axioms can be 
recast as recursion relations for certain polynomial functions which characterize the form 
factor solution completely. 

In particular, we treated the case of the free boson and the free fermion (noncritical 
Ising model with boundary magnetic fleld), where the polynomial solutions of the form 
factor axioms were shown to be identical to the explicit solutions obtained from the fleld 
theory, and it was also shown that the polynomial solutions of the form factor bootstrap 
match the full boundary operator content expected from the Lagrangian approach. 

As example for the interacting case, we flrst treated the Lee- Yang model, where the 
boundary kinematical singularity axiom makes its flrst appearance, and it is very important 
in order to distinguish between boundary conditions that have different conformal limits. 
We have also computed the spectral expansion of the two-point correlation function for the 
operator with the lowest conformal dimension and have shown that it matches perfectly 
with the ultraviolet expansion of the same correlation function obtained from boundary 
conformal fleld theory. 

Our second interacting example is the sinh-Gordon model, with Dirichlet boundary 
condition (an extension to the general case is in principle straightforward, but we decided 
to treat only the Dirichlet case to keep it short and simple). The boundary conditions with 
zero and with nonzero value of the fleld on the boundary are differentiated again by the 
boundary kinematical axiom, and we have shown that the results of the boundary form 
factor bootstrap flt perfectly well with expectations from the Lagrangian approach. 

An open question is to flnd and classify non-minimal solutions of the form factor equa- 
tions and interpret them in terms of the local boundary operator algebra of the underlying 
fleld theory, extending the method presented for the bulk sinh-Gordon model in |21| . In 
particular it is interesting to flnd out whether the counting of operators in the conformal 
limit can be matched with the full set of solutions in the interacting case. 

It is obvious that the results presented in this paper can be applied directly to any 
integrable boundary quantum fleld theory for which the factorized scattering theory is 
known, and that they formulate a well-deflned program to determine form factors and 
correlation functions of boundary operators, similar to the approach used in the bulk case. 
We have also shown how to solve the axioms for theories with diagonal bulk and boundary 
scattering. 

It is an interesting problem to extend these results to the case of nondiagonal scattering 
(with boundary sine-Gordon theory as the most prominent example). The extension of the 
axioms is straightforward: they must be decorated by multiplet indices, just like in the 
bulk case, although here we avoided to give this extension explicitly to keep the exposition 
simple. However, solving them will probably encounter much more difflculties, and just as 
in the bulk, new methods must be devised for the task, like the boundary extension of the 
Lukyanov free fleld representation in [T^ . 

The comparison to the Lagrangian and perturbed conformal fleld theory description 
would be greatly facilitated by establishing sum rules for the spectral representation of 
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the boundary correlators, similar to the c-theorem [22] and A-theorem [231 the bulk 
case, and is one of the most important problems left open by the present work. Another 
promising open direction is to consider possible applications of boundary form factors and 
correlation functions in the area of boundary quantum field theory and condensed matter. 
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A Heuristic derivation of the FF axioms 

We present some heuristic arguments, along the lines of [3], for the derivation of boundary 
form factor axioms using the boundary reduction formula 
We analyze the analyticity properties of the form factor 

as a function of the variable 6i . 

We follow the notations of jSj: The asymptotic creation/annihilation operators can be 
expressed in terms of the free asymptotic fields as 

ainiO) = 2i dxcos{k{e)x)e^^'^' dt^in{x,t) (35) 

J — oo 

a+(e) = -2i dxcos{k{9)x)e-'''^^'^'dt'^in{x,t) . 



where k{9) = msinh6' and uj{9) = m cosh 9. The in state is a free state and we have 

(010(0)1^1, 02, ... , 9n)in = {0\O{0)al{9i)\02, . . . , 0„)m (36) 
We use (jHKll together with 

OiO)^inix,t) = [O(0),<l>i„(a;,t)] + <l>,„(a;,t)O(0) 

to obtain 

= disc. - 2z / dxcosik{9,)x)e-''^^'^^'dt (0|[O(0), $,„(a;, t)] |02, . . . ,0„)m . (37) 

J — oo 

where the disconnected part is 

disc. = (0|a+ (0i)O(O)|02, . . . , 9n)in = (0 1 a+ (^i) 1 0)Ff_i (^2, ■■■,9^) 



33 



Note that in theories with nonzero vacuum expectation values of the field $ the matrix 
element (0|a^(6'i)|0) is nonzero and can be written as 

(0|a+(ei)|0) = |27r5(^^i-^) 

which corresponds to the one particle term in the boundary state in the crossed channel 
[H]. It was conjectured in and later confirmed using TBA arguments [23] that the one 
particle contribution to the boundary state has a coefficient equal to | rather than g as 
suggested in In the channel we use here this translates directly into the equation above. 

Supposing that the in field can be expressed in terms of the interacting field as t) —>■ 
Z^/'^^in{x,t) for t — s> — oo and that [C(0), $(a;, 0)] = 0, the connected part can be written 
in the form 

conn. = tZ~^/^2 j dx j rft9t|cos(A;(^i)x)e~^^(^^)*at (0| [C(0), <l>(x, t)]|^2, • • • , ^n), 

Performing the usual partial integration while taking care of the surface term we obtain 

conn. = j fc-M^i)* cos(A;(^i)a;)e(-t)(0| [0(0), J(x, t)]|^2, • • • , ^n)i„ (38) 

where J{x, t) = {D+m'^+6{x)dx}^{x, t) and the integration goes over the entire spacetime. 
The range of the integration is the interior of the past light cone due to the presence of 
9(— t) and of the vanishing of [C(0), J(x, t)] on space-like intervals. The analytic properties 
of the integral are determined by the exponent for large negative times. The exponent 
decreases if '^m{uj{9i)) > thus the in form factor {9i > 02 > ■ ■ ■ > On > 0) allows an 
analytical continuation into the domain: 

< 55m(^i) < TT 

Repeating the same procedure for the out matrix elements 



F^{-0U -02, -0n) = {0\O{0)\ - 0u -02, -0n) 



out 



we obtain the domain of analytical continuation: < '^m{—0i) < vr. 

To derive the crossing relation we consider the following matrix element 

Fin_i ■■= F^^_^{0i\02, ...,0n)= in{0l\O{O)\02, 0n)in 

Applying the reduction formula to the particle with rapidity 0i (f35|) we obtain 

= disc - 2z / rfa;cos(A;(^i)a;)e^"(^^)*9i (0|[O(0),<l>,„(x,t)]|^^2,...,^n)^„ 

J — oo 

where the disconnected part (supposing 0i > 02) is 

disc. = (0|O(0)a+ (^i)|^2, . . . , 0n)in = rn{0l\02) vnF^_2{0Z, • • • , ^n) 
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Performing a partial integration as before the result for the connected component is 

conn. = J dW"^^^)* cos(A:(^i)a;)e(-t) (0| [C(0), J{x, t)] |^2, • • • , ^n)m (39) 

which has an analytic continuation for 

-71 < Qm{ei) < 

Comparing (f38|l with (|39|) and using that mcosh(6'i+27r) = — mcosh^i the crossing relation 
^ is proved. Similar result can be obtained for an out state and the —9i < —62 < ■ • • < 
—6n < range of the parameters. 

The reflection property (Axiom II) can be shown by considering 

{01, 02,..., On) = {0\O{0)\el,e2,...,en)^n 

and crossing all particles except the one with rapidity 9n to the left. Now inserting a 
complete set of out states we have 

m(- • • |C^(0)|^„,)m = in{- ■ ■ \^iO)\fT') out out{n\ 

out 

where only the flrst two terms are nonzero: 

,„(. . . \Om9n)in = in{. • • | A(0) 1 0) (0 1 + $^(. • • \A{0)\e) out out{0\0n) ^n 

e 

The connected part gives the required R factor while the disconnected one combined with 
the disconnected part in (|H7|l and the permutation property gives the boundary kinematical 
singularity. 

The permutation property in the bulk case is usually derived from very similar argu- 
mentation we used above for showing the reflection property. Note, however, that the 
same result can be obtained from the analysis of the singularity structure of the Green 
functions: the part, which is responsible for the discontinuity in the form factor by chang- 
ing two neighboring rapidities, is related to the bulk S'-matrix. The permutation property 
in the boundary case (Axiom I) can be derived only from the second approach, namely from 
a detailed investigation of the singularity structure of the Green functions. By extending 
the result on the two point function in Pj one can show that multi point functions have 
momentum preserving parts identical to their bulk counterparts and exactly these parts 
contribute only, when two neighboring (both positive) rapidities are changed, and cause 
the same discontinuity in the form factor we met in the bulk case. 

The kinematical singularity equation (Axiom IV) can be obtained (using the permu- 
tation and reflection axioms) from the analysis of the disconnected components in the 
crossing relations as obtained for the in and for the out states: 

FZ.,{±ei\ ±02,..., ±^„) = (m ± e,, ±02, . . . , ±^„) + 27,6(01 - e2)F^_^{±es, . . . , ±^„) 

Although our derivation of the boundary form factor axioms is heuristic to some extent 
we expect that the formulation of the same ideas in a rigorously deflned quantum fleld 
theoretical framework would lead to the proper and mathematically founded derivation 
(but note that this has not been performed yet in the bulk case either). 
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B Formal derivation of the FF axioms 



Here we show how the boundary form factor axioms can be formally derived from a bound- 
ary analogue of the Faddeev-Zamolodchikov algebra. 

In the bulk case one introduces creation Z*{9) and annihilation Z{9) operators corre- 
sponding to asymptotic states. They are defined for real rapidities 9 E R and satisfy the 
following defining relations 

Z*{9i)Z*{92) = S{9^-92)Z*{92)Z*{9i) 
Z{9,)Z{92) = Si9^-92)Z{92)Zi9^) 

Z{9i)Z*{92) = S{92-9i)Z*{92)Z{9i) + 2ii5{9i- 92) (40) 

Extending Z, Z* to imaginary rapidities (treating 6' as a 27ri periodic complex variable) we 
encounter singularities in their products at 6*1 = 6*2 ± ivr with residues 

-i Res Z*{9i)Z*{92) = 1 

91=92 + 177 

and 

-i Res Z{9i)Z{92) = 1 

9i=92+in 

These can be formulated by postulating the crossing property 

Z{9) = Z*{9 + t7T) 

and taking into account the defining relations (|^ . We note that using this identification 
all of the defining relations ((i(H) can be combined into a single one 

Z*{9i)Z*{92) = S{9i - 92)Z*{92)Z*{9,) + 2n6{9, -92- in) 

In the boundary case the generators Z, Z* are defined only for positive values of the ra- 
pidity arguments and additionally two new formal generators are introduced creating the 
boundary vacuum as follows 

10)^ = 510) , b{0\ = {0\B 
We introduce two new relations 

Z*{9)B* = R{9)Z*{-9)B* 

and 

BZ{9) = BZ{-9)R{-9) 

which describe how we can extend the generators for negative rapidities. By analytically 
continuing in the rapidity again we have singularities in the operator products 

-i Res Z*(9)B* = -B* 
36 



and 

i Res BZ(9) = -B 

e=-jf 2 

They would correspond to particles with real rapidity in the crossed channel (exchanging 
time and space coordinates). These new relations can again be summarized in a single one 

z*{e)B* = R{e)z*{-e)B* + 27:6(6 - 

together with its formal conjugate 

777 n 

BZ{9) = BZ{-e)R{-e) + 2'K5{e + —y-B 

We claim that representing the form factor of the boundary operator C(0) as 

(^1, . . . , = (o|fi 0(0) z*{e,) . . . z*{e^)B*\Q) 

and supposing locality 

[0(0), Z*(^)] = 

we can recover all the non-singularity type form factor axioms immediately. For deriving 
the singularity axioms we have to observe that singularity appears not only from a single 
term. E.g. in the boundary kinematical singularity axiom, the form factor F^{6i, . . . , 6n) 
exhibits a singularity in 9i at i| coming from two places: the operator product of both B 
and B* with Z*{9i) is singular. Supposing that they appear in additive terms of the form 
factor we can obtain the desired formula. 

Finally, we note that formulating the boundary FZ algebra in the spirit of |25 might 
lead to a more rigorous derivation of our axioms. 
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